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We present an exact solution of a ID model: a particle of incident energy E 
colliding with a target which is a ID harmonic "solid slab" with atoms in its 
ground state; the Hilbert space of the target is restricted to the (A^ + 1) states 
with zero or one phonon present. For the case of a short range interaction, V{z), 
between the particle and the surface atom supporting a bound state, an explicit 
non-perturbative solution of the collision problem is presented. For finite and 
large A^, there is no true sticking but only so-called Feshbach resonances. A finite 
sticking coefficient s{E) is obtained by introducing a small phonon decay rate r] 
and letting oo. Our main interest is in the behavior of s{E) as i? — > 0. For a 
short range V{z), we find s{E) ~ i?^/^, regardless of the strength of the particle- 
phonon coupling. However, if V{z) has a Coulomb z~'^ tail, we find s{E) a, 
where < a < 1. [A fully classical calculation gives s{E) ^ 1 in both cases.] 
We conclude that the same threshold laws apply to 3D systems of neutral and 
charged particles respectively. 



PACS numbers: 68.10.Jy 03.65.-w, 03.80.+r 



The low energy behavior of the sticking coefficient s{E) of a particle striking a sur- 
face and being trapped in a surface bound state is still a matter of experimental^"^ and 
theoretical^"^ controversy. We consider here only a zero temperature target. For a particle 
treated classically which interacts with a classical elastic solid, it is known^° that the stick- 
ing probability, s{E), tends to 1, as the incident energy of the particle E ^ 0. In contrast, 
a particle whose motion is treated quantum mechanically can have a dramatically different 
threshold behavior for sticking; it is known^"^ that for sufficiently weak particle-phonon 
coupling, if perturbation theory is valid, a quantum particle experiencing a short range 
interaction will have a sticking probability vanishing as s{E) ~ E^/'^ due to "quantum 
reflection." 

A number of experiments have attempted to explore the low energy region where 
quantum effects should dictate the form of the threshold behavior of s{E). Nayak et al.^ 
found for ^He atoms striking a liquid-^He surface, s — > 0, in agreement with the quantum 
prediction. Similarly, results^ of low energy scattering of Ne from Ru(OOl) indicate the 
vanishing of s{E) as £■ — > 0. However, more recently, s{E) was determined experimentally^ 
for ultra-low energy atomic H striking a liquid- "^He surface, and did not appear to have 
quantum threshold behavior. Further threshold behavior apparently contrary to quantum 
predictions was found for positronium (Ps) on Al(lll) surfaces.^ We shall return to these 
experiments at the end. 

Theoretically the threshold behavior of s{E) for a particle coupled to a solid by a 
short range interaction was formally studied by Brenig^ who expressed many-body effects 
in a non-local, complex, energy dependent potential Uefi{r,r'; E). Assuming that Ues is 
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also short range in r and r' and has a well defined finite limit as — > 0, he found s(0) = 0, 
in keeping with quantum perturbation theory. The same conclusion^ ^ was reached for the 
important case of van der Waals interactions (~ for large z). 

Polarization efi'ects due to virtual (non- resonant) phonon excitations that are quanti- 
tatively important are neglected in Refs. 6 and 11. Knowles and Suhl-*^^ have shown that 
surface polarization effects increase s{E) at low energies. As a result of the polarization, 
the penetration of the particle's effective wavefunction into the surface region is increased. 
This effect is known to be in competition with quantum reflection in the determination of 
s{E) as E ^0. 

For the case of a charged particle which experiences a long range image potential, 
Martin et al7 calculated s{E) using perturbation theory and concluded (mistakenly) that 
s{E) oc E^/^ for small E, so that s(0) = 0. Their subsequent numerical calculations 
using the time dependent Hartree approximation indicated to them that s(0) 7^ if the 
particle- phonon coupling A exceeded a critical value Ac. They concluded that when A > Ac, 
polarization effects dominate the quantum reflection. 

We shall present results of an exactly solvable one dimensional (ID) model for quantum 
sticking. The model is sketched in Fig. la. It has an external particle interacting with 
a ID "solid slab." All motions are constrained to ID. The surface atom and the particle 
interact by a short range potential whose generic form is sketched in Fig. lb. We take for 
the Hamiltonian of the system^^ 



H = Hph + 'Hp + Hi, 
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(1) 



where 

Hph = hVLgO) gag, 

<? 

«'' = I^ + ''W, (2) 

Tip is the Hamiltonian for the particle moving in the static potential, V{z). Tiph is the 
Hamiltonian for the phonons in the solid; 7ii contains the particle-phonon coupling; m is 
the particle mass; D,g is the frequency of the phonon with wavenumber q; and a\ and ag 
are phonon creation and annihilation operators respectively. N is the number of atoms in 
the solid, M is the mass of a lattice atom, and a is the equilibrium lattice spacing. 

The phonon Hamiltonian describes a chain of N atoms coupled harmonically to their 
nearest neighbors and to fixed sites as illustrated in Fig. 1. (The coupling to fixed sites, 
which introduces a lower cut-off frequency Q,c, is needed in this ID model to prevent a 
well-known infrared divergence which does not exist in higher dimensions.) The interaction 
potential, V, is sufficiently deep so as to support a bound state with energy —Eh. 

We next restrict the Hilbert space to states with or 1 phonon present"^^. However 
we allow a particle-lattice interaction of arbitrary strength. 

We expand the wavefunction using 0— and 1— phonon eigenstates, V'o ^ind {V'i}" 

N 

*(^1,^2, ■■•,Zn,z) = (l)o{z) 'l/jo{zi,Z2, • ■ -jZn) + ^ (l>i{z)'lpi{zi, Z2, • ■ -jZn), (3) 

i=l 

where (2:1, 2:2, • • • , zn) are the positions of the chain atoms, i labels its modes, and z is the 
position of the particle. The following coupled system of equations results: 

N 

{Hp - E)Mz) + Yl ^o^iz) <P,iz) = 0, (4a) 
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{Hp + hQi - E) ^i{z) + Vio{z) 0o(^) = 0, z = 1, • • • , TV (46) 



where 



M^) = M^) = -^ y^n^)cos(f), (5) 

is the wavefunction for the particle in the i*'* channel. UVti is the excitation energy, qi 
is the wavenumber of the i^^ mode, and E is the total energy of the system. 

The coupled system of Eq. (4) may be reduced to an effective one-particle equation 
by first solving Eqs. (4 b) for the in terms of (j)Q. Substitution into Eq. (4 a) gives a 
single equation for the particle wavefunction in the elastic channel, 0o- 



(6) 



where 



AT 



U,s{z, z'; h) = Y, Uoiiz) Uoi{z') Gi{z, z') (7) 



and 



Ui-k^ - efe Jo ijJi-k^ -\- k'^ - %r] 
k = \/2mE/h, LOi = ^hfli, and = ^^fi . $(^, k) and ^biz) are normalized continuum 
and bound state eigenfunctions respectively for the static potential U{z). 77 — > 0+, so as 
to satisfy the outgoing boundary conditions for the scattered particle. 

The boundary/asymptotic conditions on 0o cire 

0o(-oo) = O (9a) 
= e-'^"" - R{k)e'^' , (96) 
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where R{k) is the elastic reflection coefficient. Eqs. (5) and (8) can be analyticaUy solved 
for s{E). Complete details are given elsewhere. We find that for finite N, s{E) = for 
all E. Only in the limit of AT — > oo does sticking occur, and s{E) is given as 

1 \P{k)\'\lmI{k)\ 
'^^^ - k \l-W{k)I{k)\^ 



where 



W{k) = 2 (^) J dz'dz"M^') U'iz') g{z',z";k) U'{z") H{z"), (11) 



p(a;)cos2 



T{k) = / du; J- ^ (12) 



and 



P{k) = -2ze^^W y ^ / dz'<^b{z') U'{z')xo{z';ky, (13) 

p{u) is the density of vibrational states per atom and ^ is a Green's function for the 
particle in the potential, U + Upoi, consisting of the static part and the part due to virtual 
excitations of the particle-phonon system, exclusive of virtual particle sticking; xo is a 
solution inU + Upoi, subject to the boundary condition that it vanish as 2; — > —00 and the 
asymptotic condition that it approach sm{kz + S) as 2; — > 00; S is the phase shift resulting 
from U + Upoi- 

In the low energy regime, for a finite range U, Xo{z, k) — > kf{z). The factor k is the 
manifestation of quantum refiection. The result is that regardless of the coupling strength, 

s{E)r. E'/' (14) 

For a neutral particle which asymptotically experiences a z~^ potential, the above result 
is still valid^^. 
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A charged particle asymptotically experiences a potential. Here (unlike for z~^) 
the WKB approximation for xo is valid for all k beyond a fixed z, and the particle experi- 
ences no quantum reflection. In the low energy regime, Xo{z, k) — >\/k g{z). We define as 

fc— >o 

Sn{E) the sticking probability for the bound state n. In lowest order perturbation theory 
in the particle-phonon coupling, we find 

Sn(^)^^ «n (15) 

where < < 1- 

We must concern ourselves with the convergence of the infinite summation resulting 
from sticking contributions from the infinite number of Coulomb bound states. The ampli- 
tude of high lying bound states near the surface behaves as n~^/^ as for pure Coulomb wave 
functions^. Thus the square of the matrix elements decrease as n~^, insuring convergence 
of the summation. In fact, most sticking occurs in the lowest bound state. 

In relating quantum sticking to classical sticking we want to point out two quite 
distinct quantum effects. 

1. A Debye- Waller like effect: In quantum mechanics there is a finite probability (even 
as N^oo) that no lattice vibrations are excited and hence the particle is refiected. 
Thus under all circumstances s{E) < 1. By contrast, classically, in the case of an 
attractive particle-target interaction, a finite amount of impact energy is delivered to 
the target, even when E—>0, because of the particle's acceleration by the interaction 
potential. When N^oo, some of this energy disappears to z — —oo. Thus for E 
sufficiently small, E < Emin, the particle cannot escape and s{E) = 1. 
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2. Quantum reflection: We consider first the particle striking a rigid target in the classical 
regime. The particle coming in with a low velocity, —Voo, spends a time of the order 
^res~^ in the interaction region, where zq is the range of interaction and is a mean 
speed in the interaction region. As Vqo— >0, v approaches a finite limit and the ratio of 
the time spent by the particle in the interaction region to the time spent in a spatial 
interval zq in the asymptotic region is 



Now consider the problem quantum mechanically for small incident energy. In the rigid 
target potential, assumed sufficiently short range, the particle is described by a standing 
wave, Xo(-2, k), with the properties 



where jV is an (irrelevant) normalization factor, S' is a phase shift, and f{z) becomes 
independent of k for small k and is of order 1. (This is well-known from so-called effective 
range theory-^^ and can easily be checked for a square well interaction potential backed by 
an infinite wall.) Thus the ratio of the probability of finding the particle in the interaction 
region to the probability of finding it in an asymptotic interval of length zq is 




(16) 




(17a) 



Xo{z,k) ^ J\f{kzo) f{z) 

2— >0 



(176) 




(18) 



Note the power of compared to the classical result, i?^/^: as E^O the quantum particle 



spends less time in the interaction region than the classical particle, by a power E^/'^. This 
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is the so-called quantum reflection. We can also note that for unit incident current, the 
probability of a quantum particle being in the interaction region is r^E^I'^. This is the 
physical origin of the sticking threshold behavior of Eq. (14). For a charged particle, on 
the other hand, one finds Pi/Poo has the form of Eq. (16), as in the classical case; i.e., 
there is no quantum reflection. 

Very recent experiments with low energy neutral particles striking low temperature 
targets have failed to flnd the threshold behavior of Eq. (14) which we expect because of 
quantum reflection. The sticking of neutral H on liquid '^He fllm^ was analyzed by Carraro 
and Cole^ using Ref. 11. They conclude that the energy (~ 10~® eV) was still too high to 
observe quantum reflection in this system. 

An ingenious experiment by Mills et al.'^ of desorption of low energy Ps from Al, led 
to the estimate s(0) a; 1 by means of a detailed balance argument. Because of the low 
mass of Ps, quantum reflection effects are expected to be much more pronounced than 
for incident atoms and molecules. According to Ref. 8, they would normally be expected 
at incident energies below 2 eV. The experiments explore energies down to 5 x 10~^ eV 
without signs of quantum reflection. 

Martin et al.^ offer as explanation their previous conclusion from numerical work 
that sufficiently strong inelastic coupling eliminates quantum refiection. Our model shows 
quantum reflection for neutral particles regardless of the strength of the coupling. We also 
want to point out that the numerical results of Ref. 8 (Fig. 2) for a charged particle show 
that s(0) oc for small coupling constant A, in line with our own results, and that, when 
plotted against A rather than log A, there is in our opinion no sign of an abrupt transition. 
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We therefore consider the apparent absence of quantum reflection observed by MiUs et al. 
as stiU unexplained. 

We acknowledge helpful conversations with M.E. Flatte, A. Kvinsinskii, Th. Martin, 
and H. Metiu. This work was supported by NSF grant DMR87-03434 and ONR grant 
N00014-89-J-1530. 
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Figure 1. (a) Schematic view of a particle with mass m impinging upon a ID sohd consisting 
of atoms with mass M. Lattice atoms are coupled to nearest neighbors and to fixed lattice 
sites, (b) Particle interacts with the end atom via a finite-range surface potential. 
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